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1. Introduction 

The implicit differential equations of the form 

A{x,t)x = b{x,t), X G i?", teR^ (1) 

are considered in this paper where A{x, t) is a ( n x n)-matrix and b{x, t) is a ( n 
X l)-vector function. The solution of equation (|1|) is called a singular solution if it 
includes the points (called singular ones) xo,to where the main determinant of matrix 
A vanishes. Such singular solutions pass through the regions where the uniqueness and 
existence theorem is violated. We will study behaviour of the solutions in these regions 
by the asymptotic methods. The initial condition for equation (|l]) is set to be 

x{to) = Xq, detv4(xo, to) = 0. (2) 

The surface detA{x, t) = could be called the singular surface in the phase space. 

In the special cases the system (|l]) could be a system of Euler-Lagrange equations 
defined by the Lagrangian L{q,q,t), where q,q G -R" are vectors of coordinates and 
speeds, respectively. One can denote the coupling {q, q}'^ by x G R^"' (a sign (^) is a 
transposition) and present the Euler-Lagrange equations 

d_dL _dL_^ 
dt dq dq 

in the form (1), where the (2n x 2n)-matrix A{x^t) contains the block I (identity (n x 
n)-matrix) and the block 

^ = {^},z,j = l,2,...n. 
dq^ dqidqj 

Therefore, the matrix A has the form 

'10 \ 



A{x,t) 



P 



The vector function of the right side is defined as 

Kx,t) = ( dL_^_ 9!L^ ) • 
\ dq dqdt dqdq " / 

In such case the singular surface is determined by the equation 
if = 0, 

where Hessian H is H = det A{x,t) = detd^ / dq^L. 

The singular solutions of the Euler-Lagrange equations, containing the singular 
points, appear in some physical applications. For example, as it has been obtained in the 
model of post-Gallelelian approximation of relativistic quasiclassical particle dynamics. 
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the Lagrangian is singular on some surfaces of the phase space |^. Another example 
of the singular Lagrangian appears in the model of higher curvature string gravity 1^, ^. 

The structure of the paper is the following. In the section 2 the few typical 
singularities for the case of degenerated matrix A by the multiplicity 1 are investigated. 
Section 3 deals with the simplest case for a degeneration by two. Prime mathematical 
demonstration examples are presented in section 4. In which it is present also the real 
physical example of typical singularities, appeared in the black-hole type solutions for 
the higher curvature string gravity model. 

2. Several types of the simplest singularities in the case of the 
single-valued degeneration 

2.1. Change of coordinates 

Without loss of generality one can consider the system (|ll-@) in the following form 

f A{x,t)x = B{x,t), 

\ x{0) = 0, detA(0,0) = 0, ^ ^ 

where x e R"", t e R\ A{x,t) = {aij{x,t) G ^^[i?" x R^]}, B{x,t) = {bi{x,t) G 
C^f-R" X R^]}, i,j = 1, n. Let the matrix A= A{0, 0) be degenerated by the multiphcity 
1, i.e. 

o 

rang A= n — 1. (4) 
The matrix A{x, t) can be rewritten as 

n 

A{x, t)=A+Y^ XiAi{x, t) + tAt{x, t), 

i=l 

where Ai, At are the ( n x n)-matrices. Let some elements of the matrices Ai{0,0) and 

o 

^4^(0, 0) be nonzero. We denote by Cq the eigenvector of matrix A 

o 

A €0 = 0, ||eo|| = 1, 



and by Cj (z = l,n — 1) the other orthogonal vectors. These vectors form the basis in 
space, and the following is correct 



Aci j^O, i = l,n- 1, 



<ei,ej>=6ij, i,j = 0,n-l 

(sign < -, ■ > denotes a scalar product). The coordinate vector x can be decomposed in 

Ci as 

n—1 

x = yeo+Y^ ZiCi. (5) 

i=l 



A= I 1 , (6) 
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o 

In this basis the matrix A, considered as algebraic operator, takes the form 

[0]r^ In-l 

where [0]^_i is the (1 x n-l)-matrix of the null elements , [0]"~^ is the (n — 1 x l)-matrix 

o o 

and (n-1 x n-l)-matrix An-\ is the block of A where 

o 

det An-11^ 0. 

The vector function B{x, t) reduces to a new form, defined by the basis e^, as 

B{x,t)^{C,df, (7) 



where C =< eQ,B{x,t) >,d = {di, i = l,n—l}, d^ =< ei,B{x,t) >. Taking into 
account the doubly continuous differentiability of matrix A elements, one can present 

n-l 

A =A +yA + J2 + tAt + [O^iy, z, t)]^, 

i=l 

where A,Ai,At are the constant ( n x n)-matrices and the quantity can be 
understood with preceding notation as 

m 

1=1 

m 



11^(0)11^0, ^,(OeC^[i?™], ^,j = l,m. 

O^(^), O^,... can be defined in the same way . 

We denote D(y, z, t) — detA(x, t) and introduce the new variable y by 

y^D{y,z,t). (8) 

2.2. Asymptotic form of equations 

For the first considered case it is possible to assume that 
dD 

\y,.,t=o^ r ^ 0. (9) 

Hence 

y = ry+ < D^, z > +D^,t + 0\y, z, t), (10) 



where = dD/dz{0, 0, 0) is the constant row vector, = dD/dt{0, 0, 0). In terms of 
y, variable y takes the form 

y = ky-< Dl z > -D^t) + 0\y, z, t) (11) 
r 

The matrix A{x, t) can be expressed in terms of y, z, t as 

t)=( y/det An-1 + 0\y, z, t) [0\y, z, t)]i_, 
' ' ' ' [0\y,z,t)]r' A^.,iy,z,t) 

O — 

where y4„_i(0,0,0) =An-i, det An-ij^ 0. Denoting by A the matrix of cofactors to 
the matrix A, one can write the inverse matrix A~^ as 

A-' = -A. 

y 

Therefore, equations (|^) can be written as 

y(^l^=AB, (12) 
So, A have the structure 

' ^ " I my, z, t)]r' yA-\, + [0\y, z, t)]lz\ ) ' ^''^ 
Taking into account equation (0), the derivative of y have the form 

dD- dD dD , 

Substituting the equations (|12]) to the expressions (|14D and taking into account equation 
(y) and the structure of the matrix A (0), we obtain 

yy = r det An-i C° + < /3, z > +7t + ©^(y, z, t) (15) 

where = C(0) in (^ and a G -R^ /? G i?""^, 7 G -R^ are the constants, obtained from 
the first derivatives of det74„_i(y, z, t), c, at = and from the vector matrixes 

[0^(y, 2;, t)]^_^, [0^(?/, 2;, t)]i~^ disposed in the matrix A. The other part of ([T^) we 
rewrite as 

= \0\y, 2, t)]riC° + + \0\y, z, t)\r\ 

and, further, in a form hke (|l^) 

yz = yiAC'+ a'Zi d") + {Az + et)C'' + [O'iy, z, t)]r\ (16) 

where d° = d{0) in (g) and the coefficients A G R"'\ A G {i?""^ x R"-^}, 6 G i?""^ 
can be calculated using the elements of matrix A. 
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2.3. Turning and intersection 

The equations (|TB|), (0), considered under the initial conditions yifi) = 0, z{0) = 0, can 
have two types of singularities. The first type we define as "the turning". 

If 7^ then equation (|T5D can be rewritten in the following simple form 

yy = Fo + 0\y,z,t), 
where Fq = r det An-i 7^ 0. In such situation equations ([l6|) have the form 

yz = yg'+[0\z,t) + 0\y,z,t)]r\ 

where = AC°+ An-i dP is the constant column vector, assumed to be nontrivial. 
That is the reason why the functions y{t), z{t) have the unique asymptotical structure 

y = ±^2^ + o{^t), z = gH+ [o{t)r,-' 

defined in the neighbourhood of the null point (Figure l(a.)). (The value o(-) may be 
interpreted in the generally accepted sense.) 

The second type of singularity for the case mentioned above appears under the 
condition C° = and is called "the double intersection". The equations (p^Sj), ( pJ] ) have 
the following general form: 

yy = ay+ <P,z> +-ft + 0^{y, z, t), 

yz = yAl\ d' + [0\y,z,t)]r'. 

Assuming that ?/ ~ + o(t^), z ^ + o{t^) at £, z/ > 0, we evidently obtain e = v = 1. 
Thus, y and z can be represented as 

y = r^t + o{t), z = it+\o{t)\l-^ 

where rj & R^, G -R"^^ are the parameters which must be calculated. It is obvious that 

but for the parameter rj the square equation is resulted in the form 
rj^ = a-T] + G, 

where G =< (3, y4„_i > +7 is the constant. Taking into account the sign of 
discriminant Dis = o? — 4(7, we declare the next possible cases: i) Dis < 0. Solutions 
are absent, n)Dis = 0. One solution is present, and iii) Dis > 0. Two solutions are 
present. The case iii) means that there are two solutions passing through a singular 
point and intersecting under the nonzero angle (Figure 1(6)). 
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2.4- More complex cases 



Further we will study the case of the singularities when the condition 
us assume that 



dD 

dy ° 



0. 



is broken. Let 



:i7) 



We split the coordinate vector z G i?"'^ to the new coordinates {y,z), where the 
dimension of new vector z is decreased by 1 and becomes equal to dim z = n — 2. 
It is assumed for coordinate y that 

The part d of the vector of the right side B{x,t) decomposes to the vector {c,d), 
where dim d = n — 2. By the analogue with (0) we can write now 



y = fy+< D", z > +D^t + 0\y, y, z, t). 
The system (Jl^) has the following form: 



1^ 



/ ^ \ / detA„_i 



y 



\ ^ J 



[O 



\n-2 



921 



\[0 



912 

yk2 + 0\x,t) 2//23 + [02ii 
yh2 + [O 



2in-2 
1 



n-2 
2in-2 



yhs + [oX-i ) 



(19) 



V d I 



where 921,912 = 0^{x,t), kj are the block components of the matrix A (|T3|) (for 
example, 1//23 is a (1 x n-2) row matrix which is a part of block yA~^^i in A disposed). 
The following relation takes place for y 



y 



dD 



dD 



oy oy oz at 



where dD/dy = 0\y,y,z,t) (see (P)) dD/dy = f + 0\y,y,z,t), dD/dz = D^ + 
[0\y, y, z, t)]i_2, dD/dt = D^ + 0\y, y, z, t). 

Expressing the variable y in terms of y by using (|^) and substituting the result to 
the system (|l^), we can write equations (O) in the general form 



(20) 



yy = C°det An~i +dy + ay+ < b, z > +btt + 0'^{y, y, z, t), 
yy = ay + ay+ < /3,z> +-ft + 0'^{y, y, z, t), 
yz = C%gy + [0\y, z, t)]?-^) + Gy + [0\y, y, z, t)]^^ 

where d e R\ a e R\ b e i?"-^, bt e R\ a e R\ a E R\ P E R''-'^, j E R\ g E R""'^, 
G E R^~^ are the coefficients calculated from the matrix A. The initial conditions for 
the system (^) are y{0) = y{0) = 0, z{0) = 0. Similarly to the situation considered 
above, the equations (^0]) have a few types of the singular solutions. Here we investigate 
some of the most realizable of them. 
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2. 5. Reflecting 



The first type appears under the condition 

CVO, 5^0, \\g\\j^O. 

We would hke to call this "the reflecting". We can flnd the coordinates y{t), 

y{t), z{t) in the following asymptotical expression in the same manner as in the above 
cases 

y^t' + o{f), y ^t" + 0(1"), z^r + o{r) (21) 

where e > 0,1/ > 0,uj > are the unknown coefficients to be determined. After the 
substitution of Eq. ( pi] ) to (|2^) and the equalization of the minimal extents, one obtains 
the following system of the algebraic equations 

e + z/ - 1 = 0, 

2^ — 1 = min{e,u,uj,l), (22) 
e + u — 1 = min{e, u, u, 1). 

The unique solution for this system is£ = 1/3, u = uj = 2/3. Let us introduce the new 
parameters 6 & R^, x ^ and A G characterizing the singular solutions by 

y = 6^+ o(^), y = xv^ + o{^), z = Av^ + [o(v^)]r^ (23) 



By substituting (|23| ) to ( pO|) and setting equal the coefficients in the terms with the 



minimal extents of t, we obtain the equations for this coefficients 
6x = 3C°det An-u 

2x' = 3a6, (24) 
2xA = SC^gS. 



System (p^ has obviously the unique solution 



X = \ -aMet An-u 5 = -V, A = g. 

V 2 3a X 

So, this type means that there is only a single solution passing through a singular 
point. The phase point moves at this solution to the singular surface y = and after 
the contact with it goes away at the same side (Figure 1(c)). 

2.6. Complicated branching 

The second case deflned by (0) appears under the condition 

CO^O, a = 0, II^II^O. (25) 



It is called "the complicated branching" . For the coefficients e, u, u we must write the 
system like ( ^2]) 



e + u-l = 0, 

2z/ — 1 = mm(2£, z/, c<j, 1), (26) 
e + uj — 1 = min{e, u, u, 1). 

We assume here that the second equation of system ( ^0]) has the form 

yy = ay+ < l3,z> +-ft + af + y- 0^{y, z, t) + 0^{y, y, z, t), 

where the coefficient a, which is assumed nonzero, can be found from A in (|T3p like 
ajjSj'j. This assumption must be added to (]T7|). By solving the system (^^, we obtain 
the unique solution e = 1/4, u = 3/4, u = 1/2. Therefore, the functions y(t),y(t) and 
z{t) have the asymptotic form 



y = 5<f7t + o(v^), y = xv^^ + o(v^=^), z = AV^ + [o(v^)]r' (27) 

where 6,x,A. are the unknown coefficients and sign s = +1 or — 1. This sign indicates 
the domain of definition of the solutions. After substituting the equations (pT]) to the 
system ([20|), like in the previous case (plf), we obtain the system 



Sx = 4sC°det An-i, 

3x^ = Mc^S^+<P,^>), (28) 
XA = 2sC^g. 

The fourth order equation for the coefficient x follows from the above system. It looks 
like 

^X^-2sC°</3, (7 > X - sd(4C°det A„_i)2 = 0. (29) 

The other coefficients 5 and A can be found consequently by the formulas 

, 4sC°detl„_i , 2sC° 
6 = , A = g. 

X X 

Now one has the question about the number of real roots of the above equation. 

For answering this question one can use the well known Shturm's method The 
algebraic form (^) is equivalent to 

Pix) = x' + aiX + ao, (30) 
where the constants ai, oq are 

ai = -^sC° < p,g >, ao = -^sa{8C^det An-iY- 

o o 

In order to find the number of real roots, we must construct the Shturm's sequence and 
then count up the difference between the number of the alternations of the signs in this 
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sequence at x = ±00. Using the Shturm's algorithm, we obtain that the number of 
roots of the equation (^) depends on sign of the expression 
R-n^ 2563 

When R < 0, roots do not exist. When R = only one root (|30|) can exist and when 
R> two roots can be found. In our case ai = sdi, Oq = sclq and sign s can be positive 
or negative. 

Therefore, we can declare the next possibilities: 

^ 4 256 , ,3 4 256 , ,3 ,4 256 , ,3 

a) a^< — \ ao\ , b) > — \ \ , c) = — \ ao \ , (31) 

for which the various types of the singular solutions of the differential equations (pO|) 
can result. The existence of the real singular solutions for the item (pl|a) determines 
the sign s as 

s = sign(a). 



In this case the singular solution of equations ( pO]) has the "turning" -like type and can 
exist only at t > or t < (Figure l{di)). 

In the item (pl|b) the sign s could be chosen positive or negative. The solution like 
the "turning" exists for each of the possible signs. Each solution contains two branches 
defined by two roots of the equation (|29|). This singular solution, placed on both sides 
of the point t = (Figure 1 (0^2) )• 

The item (0c) differs from (pT]b) with R is equal to zero at sign(t) = s = — sign(a). 
There is only one branch defined by a single root of equation (p9|) for this sign t. There 
are two roots ( |29D in the other side of t = at sign(t) = +sign(d), so, there are two 
branches of the singular solution. So, in this case three branches pass through the 
singular point. 



2.7. Triple intersection 

The third case defined by ([T7|) appears under additional condition 
C' = 0. 

It can be called "the triple intersection". The differential equations (|20|) are now 
rewritten in the form 

{yy = ay + ay+ < b, z > +btt + 0^{y, y, z, t), 
yy = ay + ay+ <P,z>+-ft + 0\y, y, z, t), (32) 
yz = Gy+[0\y,y,z,t)]r'. 
It is easy to establish that e = u = u = 1 using the analysis of extents. Therefore, we 
have the right to suggest the following behaviour 

y = St + o{t), y = xt + o{t), z = At+ [o(t)]r\ (33) 
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where obviously 

A = G. 

The unknown coefficients 6, x can be obtained from the algebraic equations 
r = aS + ax + F, 
\ x'^ = a6 + ax + f, 

where F =< b,G > +bt, f =< P,G > +7. We do not pay attention to any degenerate 
case of solution (0). In general, the system (|3^) is equivalent the algebraic third order 
equation. Hence, it can have three (Figure 1(e)), two or one root. 



(34) 



3. The double-valued degeneration 

3.1. Equations in neighborhood of a singular point 
The double-valued degeneration means that 

o 

rang A= n — 2. (35) 
Most generally in this case, there are two independent egenvectors 61,62 corresponding 

o 

to the zero egenvalue of the matrix A. Let us denote ei{i = 3,n) other basis vectors 
orthogonal to 61,62. As in (|^) 



X 



em + 62^/2 + Yu^i 



i=3 



U 

The (n x n) matrix A has the representation 



[0]r' An-2 



(36) 



n-2 



(37) 



where det 0. So, the matrix A[x,t) has the block form 

/ All A12 r^l/ ,x-|2 

A( ,^ I \ \ [0\yuy2,z,ty^ 
A[x,t) = A21 A22 

V [0\y^,y2,z,t)]r' An-2 +[0\y^,y2, z^t^nZl J 

where Xij = 0^{yi,y2, z,t). The vector of right part B{x,t) reduces to the form 
{Ci{x,t), C2(x, t), (i(a;, t))^, where Ci,di =< ei,B >, i = l,n. The main determinant 
detA{x, t) is 

deiA{x,t) = \ - det An-2 +0^{x,t), 

where A = A11A22 — A21A12. Let us choose the orthogonal egenvectors 61,62 in order to 
present the following form for the value A 

2 

A = Wi + A22/2 + E y^C>\z, t) + 0\z, t) + 03(yi, 2/2, t) (38) 

i=l 
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(it 
A is 



always could be performed by the rotation of ei, 62). The matrix of the cofactors to 
s 

(39) 



det An-2 +[0^]l [0\yu I/2, z, t)]l_^ ^ 



^ ( A22 —-^12 

A{x,t) = V ^-^21 All 

and the differential equations ([T^ ) are rewritten as 



A 





yi 


\ 


( 




m 




= A 


V 


z 


) 


V 



C2iyi,y2,z,t) 



So, 



where 



' Xyi = (A22C? - Ai2C2°)det An^2 +0\yuy2, z, t), 
= (-A21C? + AiiCO)det An-i +0\yi, 2/2, z, t), 

Xz =[OHy^,y2,z,t)Yr'■(^^o^ 

+ [0\y^,y2,z,t)]r^ 

Cf = Ci(0, ...), CI = C2(0, ...), = d{0, ...). 



+ A An-2 d° 



(40) 



3.2. The simplest case 

We'll discuss here only one case appearing at 

I I + I C° 1^ 0, (41) 

and call it "the linking of turnings" . 

It is easy to see that the variables yi,y2,z have the following unique structure 

y^ ~ t'/', Z^t. 
By defining the system (^D|) more concretely, we write it as 

Ayi = aiyi + biy2 + 0^{z, t) + 0^{yi, y2, z, i), 

Ay2 =a2yi + h2y2 + 0^{z,t) + 0'^{yi,y2,z,i), 

0-1 (42) 

Ai = giiyl + gi2yiy2 + 922yl + A a„_2 d^ 

+ eLi y^[o\z, m-' + [02(z, t)]--^ + [03(yi, ^2, ^, ^)]^^ 

where Oj, 6j G i?^, (^jj G are the coefficients which can be obtained from actual form 
of Ajj, from the — blocks of the matrix A (^) and from C°, i,j = 1, 2. Thus the 
following expressions take place: 

yi = r]iV^ + o{V^), y2 = V2V^ + o{V^), z = At + [o{t)]^-^ (43) 



13 

where the coefficients rji G R^,A G must be calculated. The coefficient s is the 

sign and can be equal only +1 or -1 like the case (|T7|). By substituting (^) to ( ^21) and 
setting equal the coefficients at the terms with minimal extents of variable t, we obtain 
the algebraic system 

f {\ir]l + \2r]2)Vi = 2s(air7i + bir]2), 
\ {Xiijj + \2ril)r]2 = '2s{a2r]i + 62^72) • 

If the values of rii,ri2 was known, the coefficients A can be determined by the formula 

^0 , i9nVl+_9i2Vir]2±922V2 
(Air/f + X2i]i 

followed from the third equation (|42|) and (^). Let assume that bi is not equal to zero 
(or 02 7^ 0). The relation between rji and rj2 has a form of square equation, as follows 
from equations (|^): 

hvl + («i - ^2)^1^72 - a,2r]l = 0. (46) 



2 , ^ _o„.„. , (44) 



A =a:.2 + ^""'^T^f^^^f (45) 



It means that real roots of equation ( ^Bj ) can exist under the condition Dis > 0, where 
Dis = («! — 62)^ + 46ia2 is the discriminant of equation (^Bf). Thus the relation is 



V2 = g^vi, 



where = (— ai + 62 ± yDis)/ (2bi) and the coefficient 771 can be determined by 

Vl ■ (Ai + Ha^f) = s{a, + 62 ± v^). 
Let us satisfy the condition 

Ai + A2((7±)VO 



which must be added to the condition (|4T| ) with the restriction | 61 | + | ^2 It^ 0. 
In this way, sign s must be chosen as 



s = sign(i?- 

where = (ai + 62 ± VDis)/{Xi + A2(5'^)^) and by joining of them, the singular 
solution has the next kind 



y2{t) 



n-2 



z = (^A„_2 d' + (Ai + A2(^7±)^)-^ ■ (^711 + ^7^^712 + (^7^)'^722) jt + [oim 

The singular solution passing through the singular point is a linking of two turnings 
at Dis > 0. Moreover, if the signs of the value are identical, the both turnings 
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are placed on the same side of point t = (Figure l(/i)). Else they are placed on the 
various sides (Figure l(/2))- Only one turning exists at Dis = 0. 

We have considered here only the typical case taking place under the condition (|35|), 
but there are a few degenerate cases, for example, one of them can appear at = 
or R~ = 0. We have restricted oneself only to the general situation. 



4. Examples and discussion 
4.1. A few prime examples 

By using the prime model examples we would like to demonstrate several types of the 
singularities discussed above. For singular differential equation 

XX = t (47) 

considered at the initial condition placed on the singular surface x = 

x{to) = 0, 

we easily determine that under the condition to 7^ "the turning" takes place: 



x{t) = ±^2{t-to)to + o{^{t-to)to). 

At to = the "intersection" takes place x = rjt + o{t), where rj = ±1. It is easy to solve 
understand Eq. (|i7| ) directly. This solution is x^ = t"^ + t^. 
For the equation 

XX = —t, 

under the initial condition mentioned above, we have "the turning" at to 7^ 0. However, 
at to = we do not obtain "the intersections" case, because the real roots of square 
equation for the parameter 77 do not exist. It is easy to see from the direct solution 
+ 1^ = tg, that it does not pass through the singular point x(0) = 0. 
To demonstrate so-called "the reflection" case, we'll turn into the next primary 
example of system of two differential equations 

I. (48) 
{y-x)x + y = 0, 

with the initial conditions x(0) = y{0) = placed on the singular surface D = det A = 

o ^ 

y = 0. The matrix A of system (|48|) is 
o f 0\ 

•4= 1 



and egenvector is Cq = (1,0)^, A Cq = 0. The vector function of the right part is 
B = (-1, 0)^. The vector VD is VD = (dD/dx, dD/dy) = (0, 1), so, r = (eo, VD) = 0, 
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but < Co, & >7^ 0. Thus, the next structure of the singular solution may take place in 
the most common case 

a; = av^ + o(v^), y = 6v^+o(v^). (49) 



By substituting this to the equations (^Sf ), we can obtain the following algebraic 
equations for the coefficients a, b 

-ah = —1, = 26, 
3 



which have the unique roots a = — v^, b = y9/2. The direct solution of the system 
(^) gives the functions x{t), y{t) in the implicit form 



t = e-^-(i-x + y), y = t + Y- (50) 

The first equation (^) is rewritten as 
, ^ 3\x 3\x^ 

from which it is easy to see the unique structure of function x{t), and further y{t), near 
the singular point x(0) = y{0) = 0. It coincides entirely to (^91). 



4-2. Black Hole example and discussion 

The singular solutions of system (|I|) can be very diverse. We considered just the cases 
which have required the least number of additional conditions. In our opinion, the 
smaller is the number of additional conditions, the more probably the case may occurs 
in the nature. 

In the model of higher curvature string gravity the four-dimensional black hole 
solutions are generated by the low energy string effective action: 

S = J d^xL = / d^x^ Ki(-i? + 2d^4>d''^) + Xe-'^ScB], (51) 

where R is the scalar curvature; is the dilaton field; mpi is the Plank mass; A is the 
string coupling parameter. The latter describes Gauss-Bonnet (GB) contribution to the 
action dsTf ). For consideration of the static, asymptotically fiat, spherically symmetrical 
black-hole-like solutions, the most convenient choice of metric is 

2 

ds^ = Adt^ - ^dr^ - f{d9^ + sin' 9dip^), (52) 

where functions A, a, f and the dilaton function depend only on the radial coordinate 
r. In this metric the Lagrangian, transformed to the most convenient for analysis form, 
is 



L(A',/',0',A,0,a,/)=Lo + ALGB, 



(53) 
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where Lq = ml,[A' f f + A{f'f + - Ap{(j)'y]/a, Lgb = 4:e-^'^(f)'[AA'{ff-A'a^]/a^ 
and the stroke denotes d/dr. 

The main determinant of Euler-Lagrange equations followed from (^) in the 
curvature gauge /(r) = r has the structure 

Dmain = AiAX"^ A"^ + BXA + C) where 

A = (-32)e~^*a2[(40'V2 - l)a^ml, + A • Ue'^'^ A' <P'] 

B = (-32)e2'^a^[(T20'mV' + X2e-^*a^ml, - ■ Se-*"* A' <f)'] 

C = 2a^[-a^ml,r^ + A' ■ ?,2e-^*ml,{a^ + 2AV) + X^Ae-""^ A' (j)'] 

At A = (GB term is absent) the asymptotically flat solution followed from the 
Euler-Lagrange ( Einstein ) equations is well known Schwarzschild's one. This solution 
is 

2M 

A = 1 , a = 1, (j) = const. 

r 

The domain of its definition is r = (0,+oo). This solution has "the intersection" type 
singular point at the event horizon A = 0. Here the Shcwarzschild's solution intersect 
with symmetrical solution which usually has been rejected because it hasn't physical 
sense (Figure 2 curve (c)). But in a case of nonzero, positive A, under the regular event 
horizon (A(rh) = 0), the Lagrangian (p3D has the further singular point at r = rg < rh. 
This singular "turning" point forms by passing trajectory of solution from the singular 
surface AA^A^+feAA+C = 0, which is nontrivial only at A 7^ 0. The domain of definition 
of the new solution becomes (rg, +00) and the solution is subjected to a turning showed 
in Figure 2. This singularity is absent in the first order curvature gravity and it is found 
as a surprise for the classical gravity. 

The investigation, performed in the present work, may be efficient for analysis of 
the compound system of the ordinary differential equation with the implicit linear 
higher derivatives when it is hard to obtain the direct solution. The investigation 
of the equations (|1]) by the numerical methods can be made without the significant 
difficulties only on the range of invertibility of matrix A. In order to extend the solution 
through the singular point, it is necessary to know the structure of the singularity. In 
the Appendix to Ref. the numerical method of integration the equations by 
additional parameter is discussed. It allows passing through a singular "turning" point 
automatically. 
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Figure captions 



Figure 1. The typical singularities in implicit differential equations. 



Figure 2. Metric functions A and C7 versus the radial coordinate T when the event horizon value Tj^ is 

equal to 20.0 Plank unit values (P.u.v.). The curve (a) is calculated with A = 1, the curve (b) is calculated 
with A = (Schwarzschild's solution). The curve (c) is the nonphysical branch passing from point T]^. The 
curve (d) shows Cr(r) function at A = 1. The arrows pointed Tg shows the positions of T g-singularity. 



